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Abstract: Unification of the gauge couplings would imply that time variations of 
the hne structure constant are accompanied by calculable and very signihcant time 
variations in the QCD scale parameter Aqcd- Since Aqcd is the dominant factor in 
setting the hadron masses, estimates made by simple variations of the fine structure 
constant may not provide meaningful limits. There may also be related variations 
in Yukawa couplings and the electroweak scale. Implications for the 21 cm hyperfine 
transition, big bang nucleosynthesis, and the triple alpha process are discussed. We find 
that the hrst of these already provides strong constraints on the underlying theory. It is 
emphasized more generally that time (and space) variations of fundamental couplings 
and their correlations may be a significant probe of ultra-high-energy physics. 
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Slow temporal variations of coupling constants are conceivable in many theories, 
e.g., those in which their values are related to expectation values of scalar fields |]l|. 
(Naturally, spatial variations may also occur p|.) This possibility has been of consider¬ 
able interest since hrst proposed by Dirac 0 over sixty years ago. Webb et al. [|] have 
recently reported that that the hne structure constant a displays time variation. The 
group reports a change ^ = —(0.72±0.18) x 10“® over the redshift range 0.5 < ^ < 3.5, 
where Aa is a at large redshift minus a at the present day. In this note we will dis¬ 
cuss how this result, and other constraints, fit together if the coupling constants of the 
strong, weak and electromagnetic interactions are unified. 

From a purely observational point of view, all fundamental constants, from the 
cosmological constant on down, should be viewed as potentially having spatial and 
time variation of an unknown sort. More interesting than the particular variation of 
any individual couplings would be the correlations in the variations of different cou¬ 
plings. Such correlations potentially will tell us how the variations are being induced, 
and therefore give insight into the underlying theory. In short, temporal and spatial 
variations in coupling constants may be a little-recognized probe of ultra-high-energy 
physics. 

Webb et aids data is based on the absorption spectra of distant quasars, but 
several constraints exist as well at other redshift values [^. The best laboratory limit 
is 1^1 < 1.4 X 10 over a period of 140 days ||^. The best geological limit comes 
from the naturally occurring nuclear reactor at Oklo in Gabon; recent re-analyses of 
this data 0], under the assumption that only a varies, and that ^ is constant, lead to 
a limit of ^ = (—0.2 ± 0.8) x 10“^^/yr over two billion years. 

There are other ongoing or planned observations which will provide additional 
limits. Fluctuations in the cosmic microwave background resulting from changing the 
ionization history of the universe could lead P] to a measurement with experimental 
sensitivity of [-^l < 10“^ — 10“^ at 1000. Nucleosynthesis 0 places bounds on 
the variations of a with roughly the same order of magnitude, though at a much larger 
redshift, z ~ 10® — 10^®, again assuming that the only variation is in a. 

If one assumes that ^ remains approximately constant throughout the history of 
the universe, the resultant value from the Webb et al. data is ^ ~ 10“^®/yr, which 
is two orders of magnitude above the Oklo bound although compatible with the other 
limits. There is, however, no particular reason to assume constant time variation, so 
until the data suggests that ^ is constant, there is no clear contradiction between Oklo 
and the new observations. 


The time-variation of other quantities can also be measured using absorption lines. 
Cowie and Songaila constrain X = a'^gpTJie/Mp, where the proton magnetic moment 
is egp/2Mp, from the 21 cm hyperfine line in hydrogen at ^ ~ 1.8. Similarly, Potekhin 
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et al. limit Y = Mp/rUe from molecular hydrogen clouds at 2 ; = 2.81: 
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= (0.7 ± 1.1) X 10-^ 
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The main point of this paper^ is to emphasize that the various fundamental con¬ 
stants are likely to vary simultaneously and in a correlated way that depends on the 
underlying physics, and that studies of such systems as the Ohio reactor, nucleosyn¬ 
thesis, or molecular/atomic absorption lines should take this into account. We will 
illustrate this by one especially well-motivated example, i.e., by theories in which the 
gauge couplings unify at some large scale (consistent with the data [0). Then a change 
in a will be accompanied by a much larger change in the strong interaction coupling and 
mass scale. For dehniteness, we consider a grand unihed gauge theory such as SU (5) 
with a single coupling constant ac- Similar considerations would apply to many string 
theories which compactify directly to the standard model.^ Unihed theories involving 
more complicated breaking patterns, exotic matter, higher Kac-Moody levels, or large 
extra dimensions would differ in the details but not in the qualitative features. 

Unihcation presumably occurs at some large scale Mg, typically of order 3 x 10^® 
GeV. At lower scales the group breaks down to SU{3) x SU{2) x U{1) characterized 
by couplings Oj, i = 1,2, 3, which however tend to ac as the scale at which they are 
measured nears Mq. For instance, in the MSSM (Minimal Supersymmetric Standard 
Model), the couplings run so that at the electroweak scale Mz they are given by 


^{Mz) = + bite 


( 2 ) 


where to = ^ ^ ~ 5.32 and a'^ ~ 23.3. The bi are hxed by the relative particle 

multiplets; in the MSSM one has bi = (^, 1, —3) . Between Mq and Mz the a* rather 
than a are the relevant couplings. Below Mz the standard model group SU (3) x SU (2) x 
U{1) is further broken to SU{3) x U{1)elm, where U{1)elm is the electromagnetic 
gauge group. This leads to a value of the electromagnetic coupling at scale Mz of 

{Mz) = (Mz) + (Mz) = 127.9. (3) 


Now imagine that is the vacuum expectation value of some held 0(f) which is 
a slowly varying function of time. (For example, 0(f) might parametrize the change 

^Some of the issues in this paper have been considered in somewhat similar contexts, for example, 



^There is an order of magnitude discrepancy between the observed and expected unification scales 
in the simplest such theories, but that is only a 10% effect in log Mq/Mz- This could be accounted 
for by threshold or other effects, and would not much affect our conclusions. 
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in the volume of some compact extra dimensions.) This in turn will induce correlated 
variation, at any given energy scale below Mq, in the SU{3) x SU{2) x 17(1) coupling 
constants. To hnd the relationship between the couplings, we make a second assumption 
(which we will relax later) that the changes in are dominated by the changes in 
— in other words, that corrections to to and to threshold effects, etc., are relatively 
unimportant.^ We also neglect higher order corrections and assume (correctly, in the 
majority of models) that Aa{Mz)/a{Mz) — Aa/a, where a = q;(0) ~ 1/137 is the 
conventional low-energy fine structure constant. We then have 

^ (4) 


so that at scale Mz 


Aa 
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Similarly, the SU{3) group’s coupling, usually called as instead of 0 : 3 , satisfies, again 
at scale Mz, 


Aas 3 as Aa 

_ _ 

0:5 So; a 
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he., the proportional variation in the strong subgroup’s coupling is almost six times as 
great as that of a, when measured at the scale Mz, where ^ ~ 15.4. 

The physical quantities we are most directly interested in are the associated relative 
changes in the neutron and proton masses and, to the extent this is possible to estimate, 
the changes in nuclear binding energies. The nucleon masses are largely fixed by the 
QCD scale Aqcd- The expression we are most interested in is therefore the relation of 
to ^ (up and down quark masses make only a minute contribution to nucleon 
masses). Aqcd is approximately the scale at which as diverges, {Aqcd) = 0 . To see 
how changes in this scale are related to changes in a, we must study the further running 
of from the scale Mz down to /i. This change is again given by the renormalization 
group equations 


M = Gs {Mz) - 




In 


Ji 


(7) 


where ag^{Mz) = 8.33 and = —11 -|- ^np, where np is the number of quark flavors 
light compared to fi. 

•^Clearly we are also assuming that the back-reaction of a varying aa on the cosmological constant 
is not large; we will comment on this below. 
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Using the above equation, the value of /x at which ag (/x) ^ 0 is given by 


Aqcd ~ 


exp 


2tt 


as {Mz) 


(8) 


where mh,mc are the bottom and charm quark masses, whose sub leading dependence 
on ag we are ignoring. This leads immediately to the desired relation 


AAqcd 27 r _i ^ 27 t Aag 

Aqcd 9 9as as 

where is measured at Mz. Using (10), we hnally arrive at a prediction 


(9) 
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(The coefficient has a theoretical uncertainty at the level of perhaps twenty percent.) 
This suggests that, within the framework of a unified theory and assuming the nucleon 
mass scale is set by Aqcd, a variation such as the one suggested by Webb et al. i of 
— = —0.72 X 10“® should be accompanied by a shift i 


m 


AMp 


—25 X 10 clearly a 

- y 

very signihcant correction. 

Before discussing the implications, let us consider possible variations in other phys¬ 
ical constants. These are likely to occur along with those in the gauge couplings, 
although the specihc relation cannot be obtained without a more specihc theory. How¬ 
ever, it is useful to parametrize them and consider likely possibilities. 

Only dimensionless quantities such as coupling constants or ratios of masses are 
physically significant. We take the point of view that the unification scale Mq is simply 
a reference scale relative to which other masses such as Aqcd are expressed. Indeed, 


the formula (|I0D should more properly be interpreted, in this language, as a statement 
about the variation in the quantity Aqcd/M e- Later we will have to consider the 
variation in the Newton constant Gn oc M~i^ relative to Me- 

We next consider the (running) Yukawa coupling ha of fermion a, which is related 
to its mass by = haV, where v ~ 246 GeV is the electroweak scale. (For the quarks, 
the physical mass is ~ {niq + Mqcd), where Mqcd ~ Aqcd, the latter contribution 
dominating for u and d.) We expect that since the gauge coupling is varying at the 
unification scale, the same will be true, to some degree, for the Yukawa couplings. Since 
the size of this effect is model-dependent, we parametrize our ignorance by introducing 


unknown constants Aq, where 


Aha{Mc) _ Xa Aac . Aa 

'^a • 


ha{Mc) 


2 ac 
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The running of ha is calculable in a given model. For the light quarks and leptons, 
the running is dominated by the gauge contributions, yielding 


ha {Mg) 
ha {Mz) 


n 


{Mg) 
ai {Mz) 


'IbT 


( 12 ) 


where in the MSSM the anomalous dimensions are given by bu-,i = —3, -^) , bd,i = 

(^, -3, and be-,i = (^, -3, O) . This implies 


Aha{Mz) 

ha{Mz) 


^ ^ ba-^iO-itG | 
\ i / 


Aa 
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where CKj is evaluated at Mz- One hnds — ba-,iOiitG ~ 4.0 for a = n or d, and 0.70 
for a = e. That is, there is a moderate magnihcation of the quark Yukawas from the 
gauge corrections (mainly from 0 : 5 ), but a much smaller effect for the electron. One 
might expect that this effect for the quarks gets further magnihed as one runs from 
Mz down to the quark mass. This is true for the bottom and charm Yukawa couplings. 
However, for the light quarks, the running should be taken down to a scale of order 
Aqgd- When this is done, and the variation of Aqcd is accounted for, it turns out 
there is a cancellation and the effect is actually smaller. Given the inherent theoretical 
uncertainties in our computations, we have chosen to ignore this small shift, and will 
simply use the for the light quarks. 

More important are the changes in the electroweak scale v, which are tied to the 
scale of supersymmetry breaking in most supersymmetric models. For example, in 
the traditional supergravity mediated models there are various soft supersymmetry 
breaking masses (scalar masses, sfermion masses, and other bilinear and cubic scalar 
terms), which are usually assumed to have the same order of magnitude nisoft at the 
Planck scale.^ Although these soft parameters run (and may change sign for scalar 
mass-squares), they are generally of the same order of magnitude at the weak scale. 
Typically, one hnds that (and the inverse Gp^ of the Fermi constant) scales as 
mloft/C(weak, where aweak = ftti + Oi 2 i and that the corresponding W and Z masses 
scale as TUsoft- The underlying mechanism for breaking supersymmetry in the hidden 
sector and therefore generating nisoft is nnknown. We shall parametrize our ignorance 
by introducing 

An ^ ^Aa 
V a 

■^We are assuming that ft, the supersymmetric Higgs mass, is comparable to the soft parameters, 
whether it is elementary or dynamically generated. 
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Then 


(15) 
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In many snpersymmetry-breaking models, a dynamical strong-conpling scale sets 
an intermediate snpersymmetry breaking scale M/, which in turn feeds into mgoft- 
Typically this relation takes a form nisoft ~ Mf n positive. It is easy to see 

that if the supersymmetry-breaking dynamics is also unihed with the standard model, 
with coupling constant ac at M^, then for any n 


Arrisoft Aac , nisoft 
-~-In ——— ~ ob -~ ([) - 

^soft Mg c^g Oi 


(16) 


or, in short, k ~ 70. (Note that there is no corresponding reason to expect the Ag to 
be large.) However, the coupling constant of the supersymmetry-breaking physics may 
be set in other ways, in which case k is completely unconstrained. 

One can now go back and investigate the correctness of our original assumption of 
considering only the change in and not that in to in (^). It is easily shown that 
the effect of the change in to induced by An is not small when |fi)| is of order 10 or 
larger. For example, the second and third expressions in (^) acquire correction factors 
(1 — IO^k) ~ (1 — 0.025k). It is appropriate, then, to recalculate the variation of Aqcd 
to see whether we have left out a large effect. We find 

- 34(1 + . 005k)— . (17) 

Aqcd 


where the coefficient of k is uncertain at the level of twenty percent due to ambiguities in 
the one-loop renormalization prescription. Note that the effect of k — 70 is substantial 
but not overwhelming. The uncertainties in its effect are comparable to or below the 
errors in [Q, as we will see. 

We finally mention that changes in the gauge couplings may also be correlated with 
changes in the cosmological constant A. Indeed, it has recently been argued that hne- 
tuned cancellations involving radiative contributions would be upset by a time variation 
in a, with enormous effect 0. Since the smallness of A is the outstanding puzzle in 
particle physics, it is difficult to know how it is affected. Naively, quantum held theory 
and classical gravity would suggest the back-reaction due to varying couplings would 
be enormous, but we already know that effective held theory reasoning is wrong, and 
do not know why. Since there is neither theoretical nor experimental guidance on this 
point, we choose to assume the cosmological constant is not ahected; whether this is 
correct is, in our view, still an experimental question. 
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As a first application, let us consider the quantities X = a^gprue/Mp and Y = 
Mp/nif.- For the nucleon masses, one has 

Y[p A TTifi A rr 

Mn = Mnuc + mu + 2md, (18) 


where Mnuc ~ SMqcd scales like Aqcd- oiMeim is the electromagnetic contribution 
contribution to Mp, with M^im also scaling like Aqcd- The variations of rrin^d and rUe 
are given in O- rrin^d are typically estimated to be a few MeV, so that Mp^n are 
dominated by Mnuc- Hence, we expect 


AMp AAqcd 
Alp Aqcd 


(19) 


(This would have to be modified if \a or k were enormous.) The proton and neutron g 
factors gp^n are well described in the constituent quark model, where they are Clebsch- 
Gordan coefficients, so we will ignore possible variation in gp. Assuming that \u ~ 
Ad ~ Ae, and denoting their average value as A, one then has 


^ ~ (-32 A A A 0.8k)— 
X a 

AY , , Aa 

- ~ (34 — A — 0.8k) -~ 

Y Oi 


(23 A 6) X 10“^ 
-24 A 6) X 10■^ 


( 20 ) 


where the numerical values are obtained using the Webb et al. value ^ = (—0.72 A 
0.18) X 10“^ and A = k = 0. (The theoretical errors in this computation, which are 
slightly smaller than the experimental errors but could straightforwardly be reduced 
by more careful calculation, are not shown.) One sees the large magnification implied 
by coupling constant unification. The values in (|^) are to be compared with the 

given in eq. (i. obtained at redshifts within the Webb 
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observational limits [[ 
et al. range. Consistency of (pOD and (P would require a rather delicate cancellation, 
with A A 0.8k ~ 32. 

We should emphasize that one ought not to conclude from this that the result of 
m is inconsistent with unification of coupling constants. Even if A A 0.8k lies outside 
the range 25 — 35, one can conclude only that unihcation requires that the variation 
observed by Webb et al. is caused by physics at distance scales long compared with 
Mq^, and preferentially affecting electromagnetic or electroweak physics. Our aim is 
to use 0 to constrain the theory, not the reverse. 

Let US now turn to another process which also could be strongly affected by coupling 
variations: nucleosynthesis. Consider the neutron-proton ratio in the early universe, 
hxed approximately at the temperature Tp at which the the weak-interaction processes 
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that interconvert neutrons and protons freeze out by comparison to the expansion rate. 
The relative abundance is 

/Y -(Mn-Mp) 


Nr, 


~ e 


( 21 ) 


where by ([T^) the neutron-proton mass difference — Mp ~ 1.29 MeV is given by 


Mn- Mp = md- 


aM, 


elm' 


( 22 ) 


For dehniteness, we will use the estimates rud — mu = 2.05 ± 0.30 MeV, and 
aMeim = 0.76 ± 0.30 MeV. The freeze-out temperature Tp ~ 0.72 MeV is set roughly 
by equating the weak interaction rate at temperature T to the universe’s expansion 
rate. Expressing the Fermi constant Gi? in terms of the Higgs vacuum expectation 
value V gives 

,4 




2^5 'J'2 

— ~ K - 

v‘^ Mpi 

where iF is a constant. Using (p!^ and dehning 

NMpi Aa 
= P- 


M, 


PI 


M 


pi 


a 


(23) 


(24) 


(recall we are holding Mq hxed,) we have 


ATf 


Tf 


= '^-3^ 


Aq; 

a 


(25) 


The variation in r = N^/Np is given by 


Ar 

— = -A 
r 


Mn-Mp 

Tp 


( 26 ) 


In addition to the explicit a dependence, there may be a much larger variation due to 


Meimr and possibly the quark mass difference and Tp. From (^), (0), and (^5|) we 
hnd 

Ar / Ah' \ 

- 37(1 - 0.08A - O.OlK - 0.02p) - , (27) 

^ V Q) / BBN 

where the subscript BBN indicates that the quantity is to be evaluated at the time of 
nucleosynthesis. Possible changes in r can be estimated from the '^He abundance |]K 
with typical estimates implying 

Ar 

-0.1 < — < 0.02, 


(28) 


















implying 


-0.003 < (1-0.08A- 0.01k 


0.02p) ( —^ 


< 0.0005. 


(29) 


^ / BBN 

(Note that the coefficient could be negative for A + 0.8k ~ 32.) If we assume a constant 
d/a, this implies the range (—3 x 10“^^/yr, 5 x lO^^'^/yr) for (1 — 0.08A — 0.01k —0.02p)^, 
which is consistent with the result of |^. However, other time dependences (e.g., a linear 
dependence on 1 + z) could give very different results; and we already know that, for 
constant Ohio and the results of 0 are inconsistent (but see below.) Thus it is vital 
to measure the time dependence of a at moderate z with high precision. 

Leaving aside time variation for the moment, we may recall that the actual val¬ 
ues of fundamental couplings have themselves been the subject of continued interest, 
particularly in light of the fact that a universe like ours may only be possible for a 
narrow range of parameters. This subject, commonly known by the catchword of the 
“anthropic principle” [0, bears on the issue of varying couplings. If the time depen¬ 
dence of the couplings was such as to bring them outside the range of allowed values 
for certain reactions, the anthropic principle applied at an earlier era might forbid an 
otherwise acceptable time variation. 

In particular, there are contexts in nuclear physics where high sensitivity to varying 
coupling constants has been noted. One might wonder, then, whether the observed 
variation might in turn affect the observed properties of stars, or affect the natural 
reactor at Oklo sufficiently to warrant a reanalysis of results obtained from it. 

Agrawal et al. have discussed at length, in a variety of settings, the anthropic 
principle’s constraints on the value of u, the Higgs vacuum expectation value. Recent 
extensions of this analysis have focused on the much tighter constraints placed by the 
existence of the triple-alpha process ^He -|-^ Be — C* where ^‘^C* is an excited 
state of carbon, 7.6 MeV above the ground state. The existence of this resonance and 
the non-existence of a resonance below threshold in oxygen ensures that (1) carbon 
is produced in stellar interiors and (2) the carbon is not immediately converted into 
oxygen. In other words the existence of carbon in our universe is sensitively dependent 
on the location of nuclear resonances in both carbon and in oxygen. The subject 
becomes particularly interesting for our considerations because Livio et al. 0 and 
more recently Oberhummer et al. have shown how sensitively the energies of these 
resonances depends on the basic structure of the nucleon-nucleon potential. As an 
illustration, consider carbon production in stars, which has been occurring in stars for 
many billions of years: the time variation of a cannot be such as to imply a change in 
the nucleon-nucleon potential by a quantity sufficient to shut down carbon production. 
There is some question of how much of a change in the nucleon-nucleon potential is 
acceptable, but it is apparently of order 1%. 
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A shift in Aqcd, as in equation (|T0[), is by itself not very important, because the 
main features of nuclei (nucleon masses, nuclear potential depths) shift along with it. 
At leading order, Aqcd is the only scale affecting QCD physics, so that all dimensionless 
quantities are shift-invariant. However, this scaling is broken by the quark masses; and 
in stars, it is also broken by any associated changes in the stellar environment, such as 
the core temperature, which might be induced by changes in the various parameters 
(including Mpi.) In particular, unlike most hadronic masses, the pion mass does not 
scale like Aqcd- This is because of its special role as a pseudo-Goldstone boson of an 
approximate global chiral SU{2) symmetry of the strong interactions, which is explicitly 
broken by the small quark mass terms. One has ~ (m^ -|- rud) with the pion decay 

constant /tt ~ Aqcd- At some level this effect will feed into other hadronic masses, such 
as that of the p, and modify their linear scaling, and will also affect such dimensionless 
quantities as the pion-nucleon coupling. However, the dominant effect should be on the 
range of the pion exchange potential, which effectively shifts the strength by a factor 
of —A{mT,/Mhad), where Mhad is a typical hadronic scale. If we assume for illustration 
that this must be < 1%, one obtains 


A{m^/Mhad) 

1 

r\j — 

A{mq/A qcd) 

1 

r\j — 

(—34 -|- A -|- 0.8k)—^ 

(^TT / ^had 

2 

T^q/AQCD 

2 

a 


< 0 . 01 . 


(30) 


The limit on is then of order 10“^ — 10“^ over a period of five billion years, a 
constraint certainly compatible with the Webb et al. observation. (This simple esti¬ 
mate is consistent, when appropriately translated, with the more thorough calculation 
of Jeltama and Sher , who considered the changes induced in the potential by vari¬ 
ations in n, keeping as{Mz) hxed.) A full calculation of the effect of varying m.„/Mp, 
using a realistic nuclear potential model along the lines of , and taking into account 
possible variations in the stellar environment, would be useful. 

Meanwhile, the Oklo reactor data have so far been analyzed in detail only for the 
effects of a direct variation of a. Changes in the effective nuclear potential, especially 
from the change in rriT^/Mhad-, could have a significant effect on our understanding of 
its implications. In particular, the statement that the combined results of Webb et al. 
and those of Oklo are inconsistent with constant ^ assumes that only a is varying; if 
the nuclear physics is also varying, then no hrm conclusion is possible at this time. 
A reanalysis of the Oklo reactor, aimed at obtaining more accurately the constraints 
that it imposes on the space of coupling constants, would therefore be welcome. It is 
possible that this reanalysis would put strong constraints on the variation of Aqcd/v, 
as well as on a, and these would be interesting to know. 

To summarize, we have considered the possibility that the variation of coupling 
constants is caused by physics at very short distance scales, where the gauge couplings 
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may be unified. If this is the case, one should not treat the variation of the fine 
structure constant in isolation. One must simultaneously treat the variations in the 
strong coupling; the running of coupling constants suggest these are very significant. 


Within our framework a measured variation of an expression of the form a 




for 


example, is due more to changes of ^ than of a. We have also calculated the changes 
induced in Yukawa couplings by anomalous dimensions. Parametrizing the changes in 
Yukawa couplings at the unification scale by A and in the electroweak scale v by k, 
respectively, we find the Webb et al. observations are consistent with other absorption 
line results (which depend on me/Mp) only for specific and large values for A + 0.8k. 
Large effects are also possible for the abundance of primordial ^iPe, but these are at a 
much earlier period and cannot be directly related to the absorption line results without 
a detailed model of the time dependence. A reanalysis of the Oklo reactor constraints 
and a detailed study of the triple alpha process, taking into account these types of 
effects, would be extremely useful. 

More generally, the study of temporal and spatial variations in fundamental con¬ 
stants is potentially a powerful probe of fundamental physics. Both the time-variation 
of these couplings (which can be constrained by combining, e.g., nucleosynthesis, cos¬ 
mic microwave radiation, and absorption line data), and the correlations amongst any 
observed variations at a fixed time, provide strong diagnostics of the underlying theo¬ 
retical structure. It is important that thorough and dedicated efforts to constrain such 
variation be made. 

As this work was completed, a paper by Calmet and Fritzsch appeared |^| with 
considerable overlap. Though there are differences in the details, the two papers agree 
on the main point: that coupling constant unification imples that a time variation in 
a be accompanied by a much larger variation in strong interaction parameters such 
as the nucleon mass, and that limits on time variation need to be recalculated and 
re-interpreted. 
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